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1. Introduction 

The discovery of the Standard Model (SM) Higgs boson, or its equivalent in theories beyond 
the Standard Model, is one of the main goals of the experimental program of high-energy colliders. 
At the LHC the Higgs boson can be discovered over the full mass range up to mn ~ 1 TeV within 
a few years of running. At the Tevatron, the CDF and DO experiments are now sensitive to a Higgs 
signal at m H ~ 165 GeV [p. 

The dominant mechanism for SM Higgs boson production at hadron colliders is gluon-gluon 
fusion, through a heavy-quark (mainly, top-quark) loop. The dynamics of this process is controlled 
by strong interactions, and thus studies of the effect of QCD radiative corrections are necessary to 
obtain accurate theoretical predictions. 

In QCD perturbation theory the leading order (LO) contribution to the gg — > H cross section 
is proportional to ce|, as being the QCD coupling. The QCD corrections have been computed at 
next-to-leading order (NLO) [||, [§] in the heavy-top limit, and with full dependence on the masses 
of the top and bottom quarks [[|]. Next-to-next-to-leading order (NNLO) corrections have been 
obtained in the heavy-top limit [ffl]. These QCD corrections, which are dominated by radiation of 
soft and virtual gluons [g], lead to a substantial increase of the LO result. The QCD computation 
up to NNLO has been consistently improved by adding the resummation of soft-gluon logarithmic 
contributions, up to next-to-next-to-leading logarithmic (NNLL) accuracy [B]. 

Recent years have also seen a substantial progress in the computation of radiative corrections 
to more exclusive observables [|7|, [| |9[ fTC| ]. In particular, in the case of Higgs boson production, 
two independent fully exclusive NNLO computations are now available [f| S]. 

This contribution is divided in two parts. In the first part we discuss an update for the total 



cross section [11], and review the corresponding uncertainties. In the second part we consider 



the fully exclusive NNLO calculation, and report on a study Q12| ] of the impact of QCD radiative 
corrections on the Higgs boson search at the Tevatron. 



2. Total cross section 

In this Section we present an update [|TT|] of the NNLL+NNLO computation of Ref. [g]. The 
results are obtained using the MSTW2008 NNLO partons [0|. We first consider the top-quark 
contribution in the loop, and perform the calculation up NNLL+NNLO in the large-m, limit. The 
result is rescaled by the exact m t dependent Born cross section: recent work has definitely shown 
that this procedure provides an excellent approximation (to better than 1 % for m# ^ 300 GeV) of 
the exact top-quark contribution Jl4|]. We then consider the bottom-quark contribution. Since in 



this case the effective theory approach is not applicable, we follow Ref. [ |15[ ] and we include this 
contribution up to NLO only [g]. Finally, we correct the result by including the EW effects [[if]] as 
evaluated in Ref. [17]. Our central predictions (a best ) are obtained by setting the factorization {jip ) 



and renormalization (;ii#) scales equal to the Higgs boson mass. Our results for the Tevatron and 
the LHC (y/s = 14 TeV) are presented in Tables 1 and 2, respectively. Comparing to our previous 
predictions (see Tables 1 and 2 of Ref. [Q]), the cross sections change significantly. At the Tevatron 
the effect ranges from +9% for m H = 1 15 GeV to -9% for m H = 200 GeV. At the LHC the effect 
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goes from +30% for m# = 115 GeV to ±9% for = 300 GeV, the increase being mainly due to 
the new MSTW2008 PDFs. 



The calculation discussed above is now available through an online calculator []18[], that can 
be used to reproduce the results of Tables 1 and 2 or to repeat the calculation for different Higgs 
boson masses and/or collider energies. 

Our results for the Tevatron can be compared to those presented in Ref . [15], obtained using the 
same set of PDFs. This computation includes an estimate of mixed QCD-EW contributions. The 
main difference with our work arises in the calculation of the top-quark contribution to the cross 
section. In Ref. [ 15] the latter contribution is computed up to NNLO but choosing pip = piR = m#/2, 
as an attempt to mimic the effects of soft-gluon resummation beyond NNLO. The final numerical 
differences at the Tevatron turn out to be small and of the order of a few per mille at the lowest 
masses, increasing to 2.5% at the = 200 GeV. 

The NNLL+NNLO calculation discussed above could be improved in various respects. Log- 
arithmically enhanced terms beyond NNLL from the Sudakov exponent have been evaluated in 
Ref. [19, 2(J. Their effect, when combined with a full N 3 LO calculation, can lead to a reduction of 
scale uncertainties to about 5% [19]. The exact small-x behavior of the NNLO coefficient function 



is also known Q21[ , [22| ] and could be included. The numerical effect is, however, smaller than 1% 
for a light Higgs. By contrast, the uncertainty that affects the Higgs production cross section is still 
large. The uncertainty basically has two origins: the one coming from the partonic cross sections, 
and the one arising from our limited knowledge of the PDFs. 

Uncalculated higher-order QCD radiative corrections are the most important source of uncer- 
tainty on the partonic cross section, and are estimated through scale variations. The scale uncer- 
tainty of our results (see Tables 1 and 2) is about ±9 — 10% at the Tevatron and ranges from about 
±10% (m H = 1 10 GeV) to about ±7% (m H = 300 GeV) at the LHC. We note that the effect of scale 
variations in our resummed calculation is considerably reduced with respect to the corresponding 
NNLO result. The reduction is more sizeable at the Tevatron, where the resummation effect is more 
important. 

The other important source of uncertainty in the cross section is the one coming from PDFs. 
The MSTW2008 NNLO set provides 40 different grids that allow us to evaluate the experimental 
uncertainties. The outcoming uncertainties (at 90% CL) are reported in Tables 1 and 2. At the 
Tevatron the effect ranges from ±6% (m# = 115 GeV) to about ±10% (»i# = 200 GeV), while at 
the LHC it is about ±3% in the mass range we have considered. 

A related and important uncertainty is the one coming from the value of the QCD coupling. 
Higgs production through gluon fusion starts at &{tt§) and thus this uncertainty is expected to have 
a relevant role. Recently the MSTW collaboration has studied the combined effect of PDF+a s 
uncertainties [^J. We find that at the LHC the PDF+a s uncertainty is about 7% at 90 % CL 
(rriH < 300 GeV), whereas at the Tevatron it ranges from 7 to 18% (ra# < 200). In particular, for 
m H = 165 GeV, we get at the Tevatron a bes t = 0.389 pb^|(scale)^}^(a s +PDF@90%CL). 

We finally point out that, besides MSTW, we have at present only two other NNLO parton 
analyses: ABKM09 [g] and JR09VFNNLO Q. A comparison of the central values of the 
cross section shows that at the LHC ABKM09 (JR09VFNNLO) result is smaller than MSTW 
by about 7% (11%) for m H = 115 GeV and by 11% (8%) for m H = 300 GeV. At the Tevatron 
ABKM09 (JR09VFNNLO) result is smaller than MSTW by about 26% (2%) for m H = 165 GeV. 
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Table 1: Cross sections (in pb) at the Tevatron (Hf = iir — nin) with y/s 
the MSTW2008 [13] parton densities. 
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Table 2: Cross sections (in pb) at the LHC ( )if = <1r = mu) with ^Js = 14 TeV using the 
MSTW2008 /|/5|/ parton densities. 



4 



QCD effects in Higgs boson production at hadron colliders 



Massimiliano Grazzini 



Although these three NNLO sets are obtained through different approaches, the large differences 
in the corresponding results confirm that the uncertainty in the total Higgs production cross section 
is still large and, at least at the Tevatron, dominated by PDFs. 



3. QCD effects in H WW -»■ Ivlv at the Tevatron 

In the previous Section we have discussed in detail perturbative predictions for the fully in- 
clusive Higgs production cross section. Total cross sections, however, are ideal quantities: experi- 
ments have always a finite acceptance. To properly take into account the kinematical cuts applied 
in the experimental analysis, fully exclusive calculations are actually needed. LO calculations are in 
this respect straightforward: one can compute the relevant matrix element and integrate it numer- 
ically over the multiparton phase space. Beyond LO the QCD computation is affected by infrared 
singularities that prevent a straightforward implementation of numerical techniques. 

In particular, at NNLO, only few fully exclusive computations exist, due to their substantial 
conceptual and technical complications [f7[ ||, ^ [TOQ . For Higgs boson production through gluon 
fusion, two fully independent computations exist that are implemented in available numerical codes 

ill- 

The CDF and DO collaborations at the Tevatron have recently reported a combination of their 
results up to 5.4 fb _1 . According to this combination, a SM Higgs boson of masses between 
163 and 166 GeV is excluded at 95% CL. In this mass region the signal is dominated by the 
H — > WW — > Ivlv channel, which provides a small number of events over a huge background. 
After a first cut-based selection, background processes remain dominant and processing of real 
data and Monte-Carlo simulations with Artificial Neural Network (ANN) methods follows. Given 
the sensitivity of the gluon fusion cross-section to higher order effects, it is important to establish 
that the sophisticated methods used in the Tevatron analysis account for these effects within the 
estimated uncertainties. 



In Ref.[|12|] we have performed a study of the inpact of cuts on the Higgs boson signal. We 
assume a Higgs boson mass mjj = 160 GeV and apply the following cuts. We consider the H — > 
WW — > pL + pL~vv channel and require at least one lepton with pj > 20 GeV and |tj| < 1.1. The 
invariant mass of the charged leptons should be ra« > 16 GeV. Leptons should be isolated: the total 
transverse energy in a cone of radius R=0.4 should be smaller than 10% of the lepton pj. Jets are 
defined according to the kj algorithm with D = 0.4: a jet is required to have pj > 15 GeV and 

hi < 3. 

We define the variable MET* as 

f MET , <b > nil 
MET* = <^ , , (3-1) 

| MET x sin , <7i/2 ' 

where <p is the angle in the transverse plane between the missing transverse energy, MET, and the 
nearest charged lepton or jet. We require MET* > 25 GeV, which suppresses the background from 
Drell-Yan lepton pairs and removes contributions from mismeasured leptons or jets. To suppress 
the tt background, we require at most one hadronic jet. 

With the cuts discussed above the NLO K-factor is reduced from 2.41 to 2. 15 and the NNLO 
K-factor is reduced from 3.31 to 2.59 (for jj.p = }JLr = run)- In order to study the stability of pertur- 
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bative collections in the presence of these cuts we have studied a set of kinematical disuibutions 
that can be computed with our NNLO programs. We consider the transverse momenta of the lead- 
ing and trailing lepton p^ d , P^ {t , the invariant mass of the charged lepton, mu, the azimuthal 
separation of the charged leptons in the transverse plane and the missing transverse energy, MET. 

A study of these distributions up to NNLO does not show significant instabilities. We have then 
compared the results to the same distributions obtained with PYTHIA and MC@NLO, rescaled so 
as to match the total NNLO cross section, without finding significant differences. 

An essential part of the experimental studies concerns distributions of discrimination variables, 
defined via ANNs. In Ref. [gj] we have studied for the first time an ANN output distribution up to 
NNLO in perturbation theory, using as input variables the leptonic quantities defined above. The 
results confirm the agreement discussed above. 

Despite the agreement in the shape of the leptonic distributions, confirmed by the ANN analy- 
sis, a comparison of the acceptances show some discrepancy. The acceptance obtained with HER- 
WIG and MC@NLO is consistent with the acceptance from the NNLO calculation. In contrast, we 
find that the acceptance computed with PYTHIA is between 12% and 21% smaller than the NNLO 
acceptance, depending on the choice of the renormalization and factorization scales. This result 
is not significantly altered by hadronization and underlying event effects and appears instead to be 
related to the matrix element and parton shower implementation in PYTHIA itself. 

4. Summary 

We have presented updated predictions for Higgs boson production at the Tevatron and the 
LHC, and discussed their uncertainties. We have presented a study of the impact of QCD radiative 
corrections on the Higgs search in the H — > WW — > Ivlv channel at the Tevatron, based on the 
NNLO calculations of Refs. ^]. This study shows that these NNLO programs can provide an 
essential help in the validation of the results from standard Monte Carlo event generators and in the 
assessment of theoretical uncertainties. 
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